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In agreement with a recent experimental discovery by IXi et al. ] (|2009l) . we also find a 
sloshing mode in experiments on the large-scale circulation (LSC) of turbulent Rayleigh- 
Benard convection in a cylindrical sample of aspect ratio on e. Th e sloshing mode has 
the sa me frequency as the torsional oscillation discovered by iFunfschilling and Ahlers I 
( 2004 ). We show that both mod es can be described by an extension of a model developed 
previously (jBrown and Ahlers 1 .2008a, bl which consists of permitting a lateral displace- 
ment of the LSC circulation plane away from the vertical center line of the sample as well 
as a variation in displacements with height (such displacements had been excluded in the 
original model) . Pressure gradients produced by the side wall of the container on average 
center the plane of the LSC so that it prefers to reach its longest diameter. If the LSC is 
displaced away from this diameter, the walls provide a restoring force. Turbulent fluctu- 
ations drive the LSC away from the central alignment, and combined with the restoring 
force they lead to oscillations. These oscillations are advected along with the LSC. This 
model predicts the correct wavenumber and phase of the oscillations, as well as estimates 
of the frequency, amplitude, and probability distributions of the displacements. 
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1. Introduction 

Rayleigh-Benard convection (RBC) cons ists of a fluid sample heated from below and 
cooled from above (for a review, see A hlers et al. I (2009)). In turbulent convection, a 
hot (cold) thermal boundary layer at the bottom (top) becomes unstable due to buoy- 
ancy and emits hot (cold) volumes of fluid known as "plumes" which detach from the 
boundary layer and transport heat vertically. In cylindrical containers with aspect ra- 
tio F = D/L w 1 (i is the height and D is the diameter of the sample) these plumes 
contribute to the driving and are carried by a large-scale circulation (LSC). The LSC 
on average forms a loop with up-flow and down-flow on opposite sides of the sample. 
The dynamics o f the LSC include spontaneous d i ffusive azimuthal meande ring of the 
LSC orientation ( ISun et al. ll2005bl : IXi et al. IbOOd iBrown and Ahlers 'l2006a''^, and re- 
oricntations bo th by azimuthal rotations ( Cioni et al. I|l99 7: Brown and Ahlcrs1 l2006ah 
and cessations ( Brown et al. 2005t Brown and Ahlers l2006a|). All of these phenomena 



are no w understood on the basis of a simple stochastic model (jBrown and Ahlers 112007a , 
2008al) derived from the Navier-Stokes equations by retaining only physically important 
terms, by taking volume averages that lead to two coupled ordinary differential equa- 
tions, and by adding noise terms to represent the influence of the small-scale turbulent 
background fluctuations. 
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It has long been observed that there are temperature a nd velocity osc illations with fre- 
quen c ies corresponding to the turnover times of the LSCJ Heslot et aL l[l987: Sano et al. 
1989t ICastaing et al. I 19891: ICiliberto et al. I Il996t IXakeshita et al. I Il996t ICioni et al. 



1997; 'Qiu and Tong "2000'.l2001a"b':'Niemcla et a^. 200l'' Qiu and T ong "2002HQiu etTl 
2004; Funfschilling and Ahlers 2004 ; Sun et al. 2005a : Tsuii et al. . 20051. Despite these 
many measurements the spatial flow structure corresponding to the oscillations remained 
unclear for a long time. A few years ago an analysis based on shadowgraph images first 
showed that the spacial structure responsible for the oscillations has a twisting, or tor- 
sional, component around the polar axis of the cylindrical container with out-of-phase os- 
cillating azimuthal flows in the top and bottom halves of the sample (IFunfschilling and Ahlers 
[2004). The torsion al shape was conflrmed m ore recently by temperature measurements 
along the side wall ( Funfschilling et al. 2008h . Very r ecently an additi onal sloshing com- 
ponent o f the oscillat ing structure was discovered by IXi et al. I (|2009h and described in 
detail by Zhou et al. I (,2009). This mode consists of a time-periodic lateral displacement 
of the entire LSC circulation plane away from the center line that is in phase along the 
entire cylinder axis. 



The re have been previous attempts to flnd the origin of the oscillations. IVillermaux 



( 19951) presented a model that attributes the temporal oscillations to periodic emission 
of thermal plumes from the boundary layers. However, recent experiments showed that 
there a re no periodic correlations between the thermal signals of the top and bottom 
plates (|Xi et al. 1 120091) . suggesting tha t the source of the observed oscillations is to be 
found in the bulk. Resagk et al. ( 20061 ) predicted an in-phase torsional oscillation mode 
in ellipsoidal containers. The oscillations only existed when the cUipticity exceeded a 
critical value; thus this model does not provide a mechanism for any oscillations in 
cylindrical containers where the oscillations are usually observed. 

The model proposed by us, which describes the spontaneous meandering, rotations, 
and cessat ions with a pair of physical l y motiv ated coupled stochastic ordinary differential 
equations te rown and Ahlers" I l2007all2008ah . was extended recently (jBrown and Ahlers 



[2008b) with perturbative terms to describe systems with various asymmetries. In one 
case we predicted an in-phase torsional oscillation where the restoring force comes from 
the pressure gradients due to the side wall in containers with non-circular cross-sections. 
At the time we did not anticipate the existence of perturbations in which the LSC 
plane is displaced away from the central cylinder axis, and thus we did not predict the 
recently characterized sloshing mode. In the present paper we show that the model can 
be extended to describe the physical mechanism for and the structure of the observed 
oscillations of the LSC in cylindrical containers by simply relaxing the assumption of a 
centered LSC circulation plane. We further show that the torsional and sloshing modes 
can be described simultaneously as a pair of oppositely traveling waves. This behavior can 
be predicted by an extension of the model in which we remove the simplifying assumption 
of a planar (height-independent) LSC and by adding the appropriate advection term. 

In the next section we p r esent t he results of new analyses of experimental data from 
Brown and Ahlers I (|2006al l2008"al[bl ) . First, in ii l2.1[ we discuss the sinusoidal fltting- 



function (SF) analysis method used by us previously and t he "temperature-extremum- 
extraction" (TEE) method introduced bv Ixi et alA (|2009h . There we also present an 
extended SF method which includes higher harmonics that enable the detection of the 
sloshing mode and a measurement of the sloshing amplitude. In that section a compari- 
son of time series for the slosh displacement angle obtained by the extended SF method 
and the TEE method is presented. In ii l2.2l we examine the probability-distribution func- 
tions of the slosh displacement an d of the torsional displacement. In §[3] we extend the 
model presented by us previously ( Brown and Ahlers 2007aL 2008al ) and show that it 
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can explain the existence of the sloshing mode as a consequence of the restoring force 
that arises from a lateral displacement of the LSC circulation plane due to the turbu- 
lent background fluctuations and away from the sample centerline. In §2] we show that 
this model a lso leads naturally to the torsional oscillation discovered several years ago 
(jFunfschilling and Ahlers 2004 ): the origin of this mode had been unexplained before 
because the model had been restricted to a LSC with a circulation plane that contains 
the vertical sample center line. Finally, in §[5l we summarize the contents of this paper 
and discuss the broader applicability of its approach to other geometries. 



2. Experimental results for the slosh displacement 

2.1. Analysis methods and the slosh-displacement angle 

We use data obtained originally bv lBrown and Ahlers ( 2008a ) to search for the sloshing 
mode. The results presented here are for a Rayleigh number Ra = (ig/STL^ / {k,v) = 
1.1 X 10^° and Prandtl number Pr = v/k = 4.38 (/3 is the isobaric thermal expansion 
coefficient, g the acceleration of gravity, AT the applied temperature difference, L the 
sample height, v the kinematic viscosity, and k the thermal diffusivity) . The sample was 
a cylinder of aspect ratio F = D/L = 1.0 {D is the sample diameter). The data consisted 
of temperature measurements by thermistors embedded in the side wall. There were 8 
thermistors, equally spaced azimuthally, at each of three different heights: the mid-height 
which is defined to be at z = 0, height z — L/A, and height z — — L/4. 

Previously, in order to determine the orientation Oq of the LSC, we fit a sinusoidal 
function 



T{e) = To + S cos{0 - Oq) 



(2.1) 



to th e thermistor temperatures measured around the side wall at a given height (jBrown et al 



2005t I Ahlers et 011120061 : iBrown and Ahlers l[2006al : iFunfschilling et al. l[2008l ). This fit- 



ting function was useful for studying azimuthal diffusion and cessations of the LSC, as 
well as torsional displacements and oscillations where the maximum and minimum of the 
temperature profile have a phase separa tion of tt. The sloshing oscillation was discovered 
and characterized ( Xi et a l. 2009; Zho u et al. 1 12009 1 using instead the "temperature- 
extremum-extraction" (TEE) method for obtaining the orientation ah of the maximum 
and the orientation ac of the minimum of the azimuthal temperature profile. The TEE 
consists of fitting a quadratic functions to a given temperature extremum and the two 
neighboring temperatures. The extremum of the quadratic fit is taken as the orientation. 
The sinusoidal fitting function has the advantage of using all the information from the 
eight temperatures at a given height, but does not allow the phase difference between cxh 
and ac to differ from tt. Thus it can not detect the sloshing mode. The TEE uses infor- 
mation from only three thermometers at a time, but does not have the phase constraint 
and thus can reveal the slosh ing mode. We used the TE E method to re-analyze data from 
earlier experiments ( Brown and Ahlers 2006a . 2008af bl) and found that those data also 
contained information about the sloshing oscillations. We define the angles au and Uc with 
a built-in phase difference of tt so that au = ac corresponds to a centered LSC circulation 
plane and counter-clockwise displacements are positive. The slosh displacement-angle is 
then a' = {at — ac)/'2.. An example of a time series of the slosh displacement angle 
a' obtained with the TEE method is shown as open circles in Fig. [TJ Since the TEE 
method fits to only three thermistor readings whereas the sinusoidal fit uses eight, the 
TEE method is more susceptible to being biased by local temperature fluctuations when 
the signal is very noisy [in a typical case, the root-mean-square noise amplitude due 
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Figure 1. Time series of the slosh displacement angle a' — (a^ — ctc)/2 for the middle-row 
thermistors at i? = 1.1 x 10^". Solid circles: from the Fouri er-moment method . Open circles: 
from the temperature-extremum-extraction (TEE) method of I Zhou et al. I (|2009r ). 



to tur bulent fluctuations is approximately 0.28(5 at i? = 1.1 x 10 (jBrown and Ahlers 



2007lJl ]. Thus we developed a method of characterizing the side-wall temperature-profile 



based on Fourier moments. The dominant term was the cosine fitting function Eq. 12.11 
We took the Fourier series of the difference between that fit and the actual signal. While 
cosine terms are still symmetric and thus cannot contribute to a slosh displacement, the 
antisymmetric sine terms can characterize the off-center displacement as a perturbation 
of the dominant cosine term. We calculated the sinusoidal Fourier amplitudes 



A„ = [T-To-5 cos{9 - 9o)] sm[n{9 - Bo)] 



(2.2) 



The n = 1 term was found to be zero, because all it would do is shift the tempera- 
ture profile, which had already been taken care of by the fit of 9o- Higher-order terms 
showed strong fluctuations, and the power spectra of mode amplitudes corresponding to 
n = 2,3,4 are shown in Fig. [21 These higher-order moments all have the same back- 

? round, simil ar to th e backgrounds of other power spectra measured in the same system 
Funfschillin g et al. ][2008). The ji — 2 term has a peak in the power spectrum at a 
frequency of 0.02011 ± 0.00003 Hz, equal to that observed by cross-correlations of ther- 
mistors on opposite sides of the sample which was 0.02008 Hz for this i? = 1.1 x 10^° 



([Brown et al. 1 120071 ). This suggests that the cross-corr elation also measured the slosh 



frequency, as suggested originally by [Zhou et al. ( 20091 ). Since the A2 term is the only 



one which differs from a noise distribution, the temperature profile can be described by 
the cosine term plus the n = 2 term. For A2 ^ 0, the temperature profile is tilted so that 
the extrema become closer together, corresponding to a slosh displacement. The slosh 
displacement- angle a' satisfies the equation 5 sin a' = 2^42 cos(2a'). A time series of a' 
comparing this fitting method with the TEE method is sh own in Fig. [T] as solid circles. 
From this time series, the oscillations reported bv lXi et al.] ( 20091 ) can clearly be seen. As 
expected, the modified sinusoidal fitting function is found to be smoother than the TEE 
method because it makes use of all 8 thermistor measurements. All cross-correlation func- 
tions between a' measured at any two different heights show oscillations that are in phase 
with each other, consistent with a slo s hing oscillatio n. Thus our obs erved oscillations are 
in agreement with those of Xi et al. (|2009f ) and .Zhou et al. 1 (|2009l) . 
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Figure 2. The power spectra of the sinusoidal Fourier moments of the temperature profile T{9) 
for the middle-row thermistors at i? = 1.1 x 10^". Thick solid (black online) line: 2nd order 
moment {n — 2). Thin solid (red online) line: 3rd order moment {n — 3). Dotted (blue online) 
line: 4th order moment (n = 4). The power spectra were smoothed for clarity. 
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Figure 3. Probability distributions of the slosh and torsional displacement angles obtained 
from the Fourier-moment method at Ra = 1.1 x 10^". (a): The slosh mode at z = L/4 (top 
row of thermistors, open squares, blue online), 2 = (middle row of thermistors, open circles, 
purple online), and z — —L/A (bottom row of thermistors, solid circles, red online), (b): Open 
(red online) circles: Probability distribution of the slosh displacement (a')z where {a'}z is the 
displacement averaged over the three rows of thermistors at fixed time. Solid (red online) line: 
Gaussian fit to the peak of open circles. Open (blue online) squares: the probability distribution 
of the twist displacement angle 6' — [6o{L/4:)—9t}{—L/4:)]/2. Dashed (blue online) line: Gaussian 
fit to the peak of the open squares. 



2.2. Probability distributions of the slosh and the torsional displacements 

The probability distributions of a' based on the Fourier-moment method, computed 
separately at each of the three vertical positions {z = ±L/4,0) for i? = 1.1 x 10^°, are 
shown in Fig. [3^. They do not differ very much from each othe r. They are v e ry fla t and 
even slightly double-peaked, in agreement with the results of IZhou et 0.1. I (|20n9l ). We 



found that p{a') varied depending on the fitting method used (see § 12. ip . which is why 
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Figure 4. The power spectrum P(t<j) of the slosh amphtude at Ra — 1.1 x 10^". Thin solid line 
(blue online): at z = L/4. Thick solid line (purple online): z = 0. Dotted line (red online): and 
z — —L/4. Dashed line (black online): a fit of the predicted power spectrum from Eg. l4.9l to the 
background of the mid-plane data. 



we developed the Fourier-moment method to filter out much of the noise and minimize 
the effect of local fluctuations on a'. We also computed the probability distribution 
p{{a')z) where {a')z is the average of a' over the three levels at any given moment in 
time. This is given in Fig. [31d as open circles. Although on average the slosh oscillations 
are in phase at all three heights, the shape of the height-averaged probability distribution 
changes to a single-peaked distribution that is closer to Gaussian in shape. We attribute 
this to random fluctuations of the relative phases at the three levels which are largely 
uncorrelated. A Gaussian flt to the peak gives a width of CTq.' = 0.27 rad. In §[3] we will 
present a model that yields such a Gaussian p{{a')). We will also present a model that 
includes the twist oscillation in §|31 A height-averaged twist amplitude can be represented 
by [9q{L/A)—9q{—L/4:)]/2. This is plotted in Fig.[5]3 for the same data set as open squares. 
A Gaussian fit to the peak gives a width of agi — 0.18 rad. 



2.3. Power spectra 

Power spectra P{uj) of a' are shown separately for the three vertical positions in Fig. [H 
Despite the fact that p{a') is narrowest at the m id-plane, the peak of the power spec- 



trum is larger at the mid-plane in agreement with lZhou et "on iiooi), showing that the 



oscillation amplitude is larger there. Besides showing a peak which gives the oscillation 
frequency the measured power spectrum has several distinct scaling regions. It reaches 
a constant in the limit of small frequency, scales as u;^^ before the peak, and scales as 
w"^ in the large frequency limit. In §|4]we will present a model power spectrum that fits 
well to this data. 



3. Model for the sloshing oscillation 

3.1. The basic model 

The mid-plane orientat ion of the LSC in cyli ndrical containers can be described by 
the Langevin equation (jBrown and Ahlers 112007a . .2008a) 
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(3.1) 



The damping term is due to the rotational inertia ol the LSC which is represented by 
the measured instantaneous temperature amphtude S, the average amplitude 60, and the 
damping timescale Tg. The stochastic driving term fg(t) is assumed to have a Gaussian 
white noise distribution and is presumably due to turbulent fluctuations. This equation 
describes the azimuthal meandering of the LSC orientation and, with the addition of 
an equation for d, reproduc es the more dramatic events of reorientations and cessations 



(jBrown and Ahlers I l2007a[ ). The temperature amplitude 6 fluctuates around a stable 



fixed-point value Sq. The details of the potential are not necessary for the purposes of 
understanding the slo w dynamics of the LSC but become important for the discussion for 
instance of cessations (jBrown and Ahlers Il2008bl ). The methodology used to obtain these 
equations involved approximate volume-averaging of the relevant terms of the Navier- 
Stokes equations. This determined the important physical mechanisms behind the LSC 
dynamics and gave estimates of t he measurable physical parameter s which were typically 
good within about a factor of 2 ( Brown and Ahlers 2007a . 2008al lbl). 

This model was extended bv lBrown and Ahlers ( 2008bl ) to include perturbations that 
break the rotational invariance of the original model. In particular, it provided a qual- 
itative explanation for why non-circular cross-sections cause the LSC circulation plane 
to align along the longest diameter. The predicted mechanism depends on the pressure 
provided by the side wall. The primary component of this pressure exists also for the 
rotationally invariant case and provides the force to drive the LSC in a closed loop. A 
secondary component appears when the rotational invariance is broken and pushes the 
LSC towards alignment with the longest diameter. For non-circular cross-sections, the 
side-wall pressure yields a potential proportional to D{d)^'^ where D{d) is a diameter 
function corresponding to the azimuthall y varying container diamet er for an LSC that 



i y 

always crosses the center of the container [Brown and Ahlers l (|2008bh ]. This gives a pre- 



ferred orientation aligned with the longest diameter with possible oscillations driven by 
stochastic fluctuations. 



3.2. Potential due to pressure forcing from the side walls 

Here we show that the same side-wall pressure-mechanism can explain the sloshing os- 
cillations of the LSC. The result turns out to be determined by a diameter function that 
depends on the off-center displacement of the LSC circulation plane. This function is 
analogous to the case for non-circular cross-section containers. We will calculate model 
terms for containers of height and diameter both equal to L and with circular horizontal 
cross-sections for comparison with experiments, but this geometrical restriction is not 
necessary in general (but of course the geometry has to be such that a single-loop LSC is 
a good approximation to the physical system) . To allow for a slosh displacement we must 
relax the original assumption that the LSC travels through the cell center. The diameter 
as a function of slosh displacement in a cylindrical container is D{a') = Lcosa'. For 
simplicity we consider a planar LSC with an infinitesimal width. We will cover a non- 
planar LSC in §111 The model can be simply explained based on a vector force-balance 
on the plane of the LSC at the inner edge of the viscous boundary layer at the side 
wall where the vertical velocity peaks. This result can be added as a perturbation to the 
Langevin equation Eq. 13.11 The radial pressure provided by the side wall can be broken 
down into two orthogonal components: one along the plane of the LSC which provides 
the centripetal acceleration iiy = 2u^/ D[a'), and one perpendicular to the LSC plane 
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Figure 5. A horizontal cross-section showing the side wall forcing for a horizontal displacement 
of the LSC. The solid vertical line corresponds to the plane of the LSC when it is displaced 
from a central alignment by the slosh displacement angle a' . The component iiy of the normal 
force F provides the centripetal acceleration for the LSC and the component iix acts as the 
restoring force for a sloshing oscillation. The primary LSC velocity Urf, is perpendicular to this 
cross-section. 



which must then equal iix — Uy tana'. A diagram of this is shown in Fig. [5] The LSC 
velocity at the edge of the viscous boundary layer wa s given as — so that uj^ 



and independent of a' ([Brown and Ahlers Il2008bf ). Combining these terms gives the 



contribution to angular acceleration from pressure forces 

.., _ 2ux _ 2 tang^ 

The pressure of the side wall provides a restoring force for the plane of LSC if it is 
off-center. While the LSC plane had originally been assumed to be aligned with the cell 
center, the off-center alignment can be generated by the action of turbulent fluctuations. 
We can describe the dynamics by adding a'g as a perturbative term to a dynamic equation 
for a' analogous to Eq. 13. 1[ assuming the stochastic and damping terms are the same for 
this off-center perturbation of the LSC as for other modes: 

../ Sa' 2 tana' 

^4' — + fe(^)- (3-3) 

OoTg ^cos^a' 

The slosh dynamics can be expressed in terms of diffusion in an azimuthal potential given 

by 

^.--Z-;'^"'-^- (3-4) 
where D{a') — Lcosa' . This is written in terms of the diameter fun ction with the same 



form as the effect of the cross-section geometry (jBrown and Ahlers 11200 8b). In principle 
the diameter function could be calculated for slosh displacements in different cross-section 
geometries. A more detailed derivation based on volu me-averaging the Navier-S tokes 
equations throughout the sample volume was shown bv lBrown and Ahlers 1 (l2008bh and 
produces the same result. 

We have derived equations assuming the width of the LSC is infinitesimal for simplicity. 
Accounting for a finite width of the LSC is not expected to qualitatively change results; it 
may have the effect of smoothing out the potential near the minimum, and the potential 
may increase dramatically at an angle smaller than tt/2 due to the LSC being forced to 
the side of the container. 
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3.3. Oscillations 

We can obtain a linearized Langevin equation from Eq. 13.31 by assuming small a' and 
assuming ^ » which is a good approximation in experiments except during cessations 



([Brown and Ahlers Il2008al ). This equation is given by 



a' ^ ^- - cula' + f^{t) (3.5) 
The resonant frequency of the deterministic part of the equation is 

v/^^-l/(4r|) . (3.6) 



UJr 



This corresponds to the natural frequency when the damping term is small compared 
to the restoring term, a nd is of the sam e order as the frequency of the LSC turnover as 
observed in experiment (|Xi et The natural frequency scales with the turnover 



frequency as observed because the restoring force comes from side-wall pressure which 
balances with the centripetal acceleration of the LSC. 

3.4. Probability distribution 

The probability distribution of the orientation can be obtained from the steady-state 
Fokker-Planck equation in the strong-damping limit and is 

p{a')^e^p(-^) . (3.7) 



Expanding Vg up to 2nd order gives a Gaussian distribution with widt h o-g/ = ^/D^T^juJi, 



Based on experimental measurements of the parameters Tq and (jBrown and Ahlers 
2008a[ l we expect a ctq' « 0.22 rad at Ra = 1.1 x 10^" and a Rayleigh-number dependence 
proportional to i?a~°'^^. The predicted width of the d istribution is close to the value 0.27 



rad found in § 12.21 with no significant trend with Ra ( Zhou et al. Il2009h . Since there are 



oscillations, a strong-damping solution is not exact, but nonetheless gives an estimate 
for the amplitude. 

If we consider the prediction for larger a, we note the restoring force diverges at 
a' = 7r/2, which means that the amplitude of a' is strongly bounded to remain less than 
7r/2. This is physically significant because at a' = ±7r/2 the LSC would be confined to 
a vertical line at the side of the container. The typical amplitude of oscillation is on the 
same scale, so if the potential did not turn up sharply large fluctuations could lead to a 
breakup of the LSC. This specific form of breakup has not been observed. The dropoff 
in the probability distribution shown in Fig. [3] occurs at a smaller displacement angle of 
about 0.6 rad, which would be expected due to the finite width of the LSC. 



4. The torsional oscillation 

The sloshing oscillation occurs along with a torsional oscillation with out-of-phase 
azimuthal flows in the top and bottom halves of the sample of approximately the same fre- 



quen cy and amplitude |Funfschilling and Ahlers 1120041 : 1 Funf schilling 'et 0/71120081: [Zhou et al. 



20091) . These two modes seem to be sufficient to describe all of the observed temperature 



and velocity oscillations. Even observations that at first suggested different structures 
can be described as due to projections of the sloshing oscillation combined with the 
spontaneous meandering of the LSC. These include out-of-phase a nti-correla tions be 



tween temperatures at opposite sides of the sample ([Ahlers et al. 1 12006; .Brown et al. 
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(a) top middle bottom (b) 




Figure 6. (a)The torsional and sloshing oscillations at different heights and phases in the Eu- 
lerian frame. Circled dot: au where the LSC flow is upward. Circled cross: Uc where the LSC 
flow is downward. Connecting lines: the local plane of the LSC. Arrows pointing horizontally 
from the LSC plane: restoring force on the LSC plane due to side wall pressure. Arrows be- 
tween diagrams indicate the viewing order to follow the advection of upward and downward 
moving plumes. Slosh displacements are seen at 2tj0t — 0,tt and twist displacements are seen 
at 2uj^t = tt/2, 37r/2. (b) Displacement angles of hot and cold plumes in time; they can be seen 
as traveling waves moving in opposite directions. 



20071 ). and apparent sloshing parallel to the LSC plane (jZhou et Thus if we 



can describe both sloshing and torsional oscillations, we can describe all of the observed 
oscillatory dynamics. Some measurements indicated a difference of up to 20% between 



two o bserved frequencies measured by different methods depending on Ra (jBrown et al. 



l2007l) . but for now we will consider this as a perturbation and first try to explain why the 
frequencies should be the same. Both the sloshing and torsional oscillations have been 
described in the Eulerian frame where they appear to have a phase difference of 7r/2 so 
that the sloshing oscillation amplitude is maximized while the twist oscillation amplitude 
is zero and the signal at the mid-plane follows t he bottom row on t he hot side but follows 
the top row on the cold side. It was noted by I Zhou et al. I (|2009f l that this direction is 
consistent with the direction of advection through the sample with the LSC. If we instead 
consider this flow in the meridionally moving Lagrangian frame of the LSC, the motion 
is simply a horizontal oscillation of the plumes as they are advected along with the LSC 
such that the hot and cold plumes are always moving in opposing directions. The oscil- 
lation is drawn in both the Eulerian and Lagrangian frames in Fig. [6l This suggests that 
the combination of sloshing and torsional oscillations can be reinterpreted as traveling 
waves consisting of oppositely moving plumes. 

We now take the step of writing down the equations for these traveling waves. The 
nodes of the sloshing oscillation are obse rved to be at the top and bottom plates, and 
the antinode is found at the mid-plane. (IZhou et al. 2009) The node of the torsional 



oscillation is observed to be at the mid-plane (|Funfschilling et al. I l2008f) . For ease of 



comparison to side wall measurements we will assume the path of the LSC to be square 
around the edge of the side wall with a length of AL. This differs from the assumption of a 
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circular LSC shape in In §[3] where we needed to calculate a centripetal acceleration, but 
the calculations can be treated independently and in each case are self-consistent. This 
means that there are two peaks of the oscillation per cycle of the LSC, the wavenumber 
is k = tt/L, and the oscillation frequency is 2w0 if the traveling waves are moving at 
the speed of the LSC. It has been measured that temperature and velocity oscillations 
hav e a frequency Uma^,/2L where Umax is the maximum of the vertical velocity pro- 



file (|Qiu and Tone Il2001b[ ). Particle -image velocimetry measurements have shown that 



Umax overestimate s the average LSC speed by about a factor of 1.2, i.e. Umax ~ 1-2(C^) 



Sun et al. Il2005a^. The path length I of the circulation is given by / = 3.4L for large Ra 



Sun and Xia 20051 ). so the oscillation frequency is 2.Ow0, which is self-consistent with 



the picture of two wavelengths per cycle. Adjusting the path length from the simplifica- 
tion of 4L to the measured value of 3.4i would require a small adjustment of k but does 
not shift the frequency. We can now write down the empirical equations of motion. For 
the oscillating peak positions of the hot upward-traveling waves we have 

ah ^ Acos{T:z/L-2uj^t) (4.1) 
and for the cold downward traveling waves we have 

Uc = -Acos(-7rz/L - 2uj^t) . (4.2) 

The slosh displacement is 

a' = — ~ = Acoanz/L cos{2uj^t) (4-3) 
and the twist displacement is 

e' = ^ -^sin7rz/Lsin(2w^t) . (4.4) 

This form gives a sloshing oscillation amplitude of v4cos(7rz/L) which is strongest at the 
mid-plane, qualitatively similar to experimental observations (jZhou et al. ][200i). This 



gives an amplitude of Asm{TTz/ L) for the twist which corresponds to A/^/2 at heights 
±L/4 where the amplitude is measured in experiments. We measured root-mean-square 
displacement amplitudes ctq,/ = 0.27 rad and tre' — 0.18 rad at i? = 1.1 x 10^" (see ^ 12. 2p . 
which yield a ratio of 1.5 in good agreement with the predicted ratio of \/2. 

4.1. Advected oscillation model 

Since the observed oscillations can be described by a pair of oppositely traveling waves, 
we now extend the linear oscillation model Eg. 13.51 to describe these traveling waves. To 
describe the torsional oscillation with a stochastic differential equations we must now 
allow for displacements that vary with height, i. e. a non-planar LSC. With variation in 
height z the advection term moving along withthc LSC is no longer zero as was assumed 
in earlier versions of the model (jBrown and Ahlers Il2007al ). Advection along with the 
LSC is represented by the term da = —u ■ Va « ^Uda/dz where a pathlength of 4L and 
angular frequency a;^ gives U = 2lo^L/tt. The advection is upward for the hot plumes 
and downward for the cold plumes. Taking the time derivative of this term to add it to 
a 2nd order equation adds spurious stable solutions with constant d 7^ which we will 
ignore. We rewrite Eg. 13. 51 including this term. Separating equations for au and cuc yields 



(4.5) 
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and 

ac{z) = + — + fc{z,t). (4.6) 

"Tq Zj it z 

Possible steady-state solutions of the corresponding deterministic equations have the 
form 

ah — v4exp(znfco2 — ioJnt) (4-7) 

and 

ttc — v4exp(— mfco-z — iujnt + iip) . (4-8) 

where ko = tt/{2L). For the circulation to be closed around a loop so that the plume paths 
are continuous, the solutions must have integer values of n and a relative phase difference 
tjj = {n + l)7r between and ac- Note that because we defined displacements to be 
positive if they are counter-clockwise from a top view, a positive displacement for upflow 
followed around the loop becomes a negative displacement for downflow, accounting for 
the extra phase of tt. If we consider the n = 1 mode with -0 = 0, there is an in-phase 
torsional displacement and an out-of-phase slosh occurs at the top and bottom. These 
slosh and twist displacements have a relative phase of tt/2. It turns out that the path 
lengths of the LSC with the out-of-phase slosh displacement and the in-phase twist 
displacement are exactly the same as for a centered vertical LSC up to 2nd order in a, 
so there is no restoring force for the n = 1 mode. The n — 2 mode corresponds to the 
observed oscillations of Eqns. 14.31 and 14.41 There is a restoring force against this slosh 
displacement for n = 2 because the path length is shorter than a vertical centered LSC. 
More generally, the form of the restoring force is complicated and will not give a pure 
sinusoidal motion. To get a sinusoidal solution as an approximation we calculate the 
restoring displacement with a phase difference of tt rather than across a horizontal plane 
as [ah{z) — ac{~'z)]/2 = ak.{z)[l — exp(j?/;)]/2. This displacement is equal to zero for odd 
n and equal to ah{z) for even n. While all even-order modes have some restoring force, 
the lowest-order even mode will dominate because damping will reduce the contributions 
of higher order modes. This suggests that the observed mode should be n = 2 and 
k = tt/L in agreement with the empirical equations 14.31 and l4. 41 With this restoring-term 
approximation, the solutions have the power spectrum 



Piu;) = D, -ujI- niuuj^f + [uj/t^Y] ' . (4.9) 

In the limiting case of weak damping, the resonant frequency is ujr = (1 + ^/2)uj,p. 
Resonance with a global maximum in the power spectrum occurs if l/(T^a-'^) < 2/(1 -|- 
V2). B ased on independent measurements of the model parameters (|Brown and Ahlers" 



l2008al) we obtain l/{TgUJ^) — 2.2 for i? = 1.1 x 1 0^° which appears to be slig htly above the 
resonance threshold. However, it was argued bv lBrown and Ahlers I ( 2008bl ) that, because 
of the stochastic nature of the damping term S{t) / {Sqts) of Eq. 13. H the effective damping 
is much less than the fixed-point value of 1/t^, perhaps by an order of magnitude. Those 
authors found oscillations with a natural frequency weaker by about a factor of 5 than 
for the slosh, so by comparison the slosh restoring force should easily be strong enough 
to drive the predicted traveling- wave oscillations. For any value of the damping term 
that allows resonance, the global resonant frequency remains between 2w0 and \P2lij^ so 
is within 20% of the estimated measured frequency of 2.Oa;0. 
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We fit the predicted power spectrum of Eq. l4.9l to the measured power spectrum of a' at 
the mid-plane shown in Fig.[4l We fit the background excluding the range u!r/2 < uj < 2u)r 
and allow wq, t^, n, and Dg to be fit parameters. The fit gives lo^/{2ti) ~ 0.0101±0.0002, 
27rTg = 6.1 ± 0.1 s, n = 1.95 ± 0.08, and = (1.21 ± 0.01) x lO"*. These fit values of 
Tq and Dg are within 20% of the parameter values obta ined from independent methods 
based on power spectra of 9^ ( Brown and Ahle"rs1l2008ah confirming our assumption that 



the same values are appropriate to describe the dynamics of a' . The fit value based on 
the background gives a;^ = O.SOcj^ where is the measured peak frequency, which is 
consistent with our prediction for the n = 2 mode in the strong damping limit and 
with observations. In contrast, for a linear harmonic oscillator the natural frequency is 
always greater than the resonant frequency. The fact that the fit does not produce a 
peak is not a concern because the peak is expected to be enhanced by fluctuations of 5. 
The enhancement is all near the natural frequency since LOrTg — 0.9. The region with a 
scaling of uj~^ is not expected in typical linear harmonic oscillator equation and appears 
in the model as a result of the advective term. The crossover to a scaling with a fit 
value consistent with n — 2 provides support for the advective term in the model. 

The probability distributions of angular displacements p{a') and p{0') are predicted to 
be Gaussian, see Eg. 13.71 However, the displacement distribution p{a') at a single height 
was observed to be flat or even slightly double-peaked as shown by Fig. [3^. A twist 
displacement for n = 2 has a shorter path length than a centered, planar LSC and thus 
corresponds to a lower potential. Since this effect is maximized at 7r/2 rad out of phase 
with the slosh displacement, the twist may serve to further drive the oscillation to larger 
amplitudes, increasing the effective restoring force and giving a flatter or double-peaked 
probability distribution. An average of the slosh displacements over different heights 
averages over different phases of the oscillation and thus will tend to give a smoother 
distribution. The flat probability distribution cannot be seen in the twist because there 
it is only measured based on a difference in angle at different heights. In any case, the 
probability distributions based on the height-averaged displacem ents were shown to be 
in good agreement with the model. It was observed bv Funfschil ling et al. ( 2008h that 



the torsional oscillation amplitude tends to be larger when the LSC strength S is weaker. 
This is expected to be the case both because the restoring and damping frequencies are 
physically dependent on the LSC strength and instantaneously should be proportional 
to (5, which would tend to widen the probability distribution calculated from Eq. 13.71 
when 6 decreases. A simple substitution does not give an exact result but we can make 
a qualitative prediction that the amplitude of oscillations increases as 6 decreases. 

4.2. Containers with non-circular cross-section 

We can apply this model to other geometries by generalizing the diameter function 
D{9,a'). This is valid because the potential in Eq. 13.41 has the same form as that found 
for non-circular containers where the diameter function D{9) corresponds to the length 
of a chord throu gh the cente r of a horizontal cross-section as a function of orientation 



9 (|Brown and Ahlers I l2008bl ) . Since the potential scales as D{6)^^ the potential mini- 
mum corresponds to a longest diameter. In highly elliptical or rectangular cross-section 
containers the side wall shape was predicted to induce an in-phase torsional oscillation 
that corresponds to the n = 1 mode bec a,use the restoring force is proportional to 



measured relative to the longest diameter (jBrown and Ahlers Il2008b[ ). In the interest of 



studying oscillations in non-circular containers we calculated the generalized diameter 
fu nction D(9, a') for some ge ometries . 
iBrown and Ahlersi ( 2008 b^ considered containers with elliptical cross-section defined 



by D{9) = L[l -\- ecos{29)] with ellipticity e. An elliptical container can be considered a 
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perturbation on a cylindrical container where s is the small parameter. A second order 
expansion around the longest diameter gives D{9, a') — L(l — 2ed^ — a'^/2). The n = 2 
mode should be dominant for small e while the n ^ 1 mode is overdamped, but once s 
is large enough that the n = I mode is underdamped it should soon become dominant 
because it occurs at lower frequency so the peak of its power spectrum is reduced least by 
damping. The restoring forces for the n = 1 and n = 2 modes are equal for e = 1 /4 so at 
that point we would expect the n = 1 mode to dominate but because of the complicated 
effect of fluctuations of S it is difhcult to give a more precise value for a transition. 

In a rectangular container of horizontal aspect ratio A the potential becomes asymmet- 
ric around the a diagonal. A first order expansion around a diagonal gives the diameter 
function D± m [1 — A^^{\9\ + \a'\)]Dm where the ± terms correspond to the two sides 
of the longest diameter Dm- The scaling with \a'\ rather than a'^ should give a non- 
sinusoidal os cillation and change s the expected p{a') to a double-sided exponential as 
explained by Brown and Ahlersi ((2008b ') for the in-phase torsional oscillation. The slosh 
and twist displacements provide equal contributions to the restoring force so we expect 
the n = 1 mode to dominate in rectangular containers b ecause it is at a lower freq uency. 

We note that the tilt-induced oscillations studied by iBrown and Ahlers I |2008hl are 
similar to the predicted cross-section dependent oscillations because the restoring force 
is propor tional to 9', so this restoring force can drive the n = 1 mode. The oscillations 
found bv iBrown and Ahlers ] lioosb) had frequencies smaller than uj^ when induced by 
a small tilt angle and did not have an advective scaling regime in the power spectrum 
with P{uj) (X uj~^. Thus they were best described without the advective term in the 
model. This can be the case because the restoring force is relative to fixed locations in 
the container and can perturb the LSC without the requirement of producing closed 
oscillations. We presume that either type of mode can exist for cross-section dependent 
oscillations where the non-advective mode occurs for small asymmetries and the n = 1 
advective mode only occurring for larger asymmetries. 



5. Conclusions 

We used a physically motivated stochastic model to explain the observed torsional 
and sloshing oscillations. The proposed mechanism for the restoring force is the pressure 
produced by the side wall to contain the LSC which can have a component to drive the 
LSC towards a shape with a longer path length. This result came from extending an 
earlier model used to predict that the LSC would align with the longest diameter in a 
container with non-circular cross section by simply relaxing the assumption that the LSC 
travels through the cell center. By further relaxing the assumption of a planar LSC and 
adding the necessary advective term which was zero before, this model was extended to 
describe an advected oscillation which contains both the observed torsional and sloshing 
modes. Considering pressure gradients due to the side walls, advection along the direction 
of the LSC, inertial damping, and a stochastic driving force as the relevant terms, we 
can explain the structure, frequency, and amplitude of the observed oscillations. The 
apparent torsional and sloshing oscillations can now be understood as a manifestation 
of traveling waves which are advected along with the LSC. Since the restoring force for 
both oscillation modes comes from the slosh displacement, the out-of-phase torsional 
mode cannot exist without a sloshing mode. 



The agreement of the oscillation frequency with twice the cell-crossing time (|Qiu and Tone 



2001b^ was in the past consid ered as a sign of th e alternating emissions of plumes from 



the thermal boundary layers ( Villermaux 1995f l. However, it was sho wn recently tha t 



plumes are not emitted alternately from the thermal boundary layers ( Xi et al. 20091 ) . 
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We have now shown that this agreement can be explained by an advected oscillation. It is 
notable that we never had to invoke vertical heat transport, thermal boundary layers, or 
spatially or temporally periodic plume emission either as driving forces or boundary con- 
ditions. At this point we believe we can explain the origin of all of the known dynamics of 
the LSC using stochastic differential equations based on approximate volume-averages of 
the Navier-Stokes equations. The physics has included bouyancy, viscous damping, effec- 
tive damping from rotational inertia, restoring forces from pressure gradients, turbulent 
fluctuations, and various perturbations from asymmetries. Thus the LSC dynamics seem 
to be independent of the thermal boundary layers and plumes. The dramatic cessations 
were described by a different regime of the model where fluctuations in 6 are important, 
indicating that the two sets of dynamics are independent despite the fact that the dura- 
tion of a cessation is comparable to an oscillation period. The connection is that these 
timescales were both fundamentally set by advection. 

While we calculated and tested results for cells with circular cross section, this model is 
formulated in such a way that it can be applied to cells with other geometries. Specifically, 
the diameter function D{9,a') describing a horizontal LSC cross-section is used as an 
input for the model. This can in principle be applied to other container geometries as long 
as the flow is a single convection roll. Since this model has been so siicccsshil at describing 
the dynamics with a single-roll LSC, the next logical step is to attempt to apply it to more 
complicated buoyancy-driven convection problems, many of which involve multiple-roll 
LSCs. 
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